First-principles study of electron transport in few-electron open quantum dots by the 
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Electron transport properties of few-electron open quantum dots within the spin-restricted 
Hartree-Fock approximation are studied. The self-consistent numerical calculations were performed 
for a whole device, including the semi-infinite leads, without employing any phenomenological or 
adjustable parameters. Inclusion of the non-local Fock potential brings qualitatively new physics in 
comparison to the Hartree approach: electron screening decreases, resonant energy levels become 
less pinned to the Fermi energy and clearly correlate with conductance peaks. When coupling be- 
tween the dot and leads decreases the number of electrons inside the dot becomes quantized and 
the model predicts the Coulomb blockade of electron transport. This is confirmed by comparison 
with the master equation approach for an equivalent quantum dot. 
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I. INTRODUCTION 

Electron transport in mesoscopic structures is usually 
well described within the mean-field approaches that in- 
clude electron-electron interaction at various levels of 
sophisticatioujir— The simplest model accounting for 
classical Coulomb repulsion between charged particles is 
the Hartree approach. While being the simplest it is able 
to capture many transport phenomena that are inacces- 
sible for the single-particle approaches or the Thomas- 
Fermi model. For example, conductance statistics in 
open quantum dots at low temperatures was revealed 
to be governed by the pinning effect of resonant energy 
levels 1^ However many phenomena observed experimen- 
tally demand further improvements of theoretical study 
such as inclusion of spin and of the exchange interaction. 
The hallmark example is the Coulomb blockade effect in 
quantum dots.^^ It is usually described by a phenomeno- 
logical approach involving a set of adjustable parameters 
along with completely disregarding the dot geometry;^^— 
As a result it is difficult to relate the employed parameter 
sets to the physical processes they represent and to the 
device topology studied. A model that accounts for both 
the dot and lead geometry and treats the whole device on 
the same footing in the Coulomb blockade regime as well 
as accesses direct quantitative description of Coulomb 
blockade in dots of general geometry has not been accom- 
plished yet. Such a model would allow one to capture the 
microscopic physics of electron transport in the quantum 
dots. 

Attempts to address Coulomb blockade, at least at a 
model level, are made within time-dependent density- 
functional theory (DFT)i^ii^ Electron transport is con- 
sidered to be intrinsically non-equilibrium that evolve 
in time and Coulomb blockade occurs in the long- 
time limit as a periodic sequence of charging and dis- 
charging events. Electron tunneling through the in- 
teracting many-body system can also be considered 
within GW-approximation in the framework of non- 
equilibrium Green's functions (NEGF).^-'^ By introduc- 



ing quasi-particle electronic structure it brings substan- 
tial improvements in comparison to standard NEGF and 
agreement with experiment becomes betterfli However, 
whether time-dependence is necessary for steady-state 
regime and whether electronic correlations along are re- 
sponsible for Coulomb blockade are not clear. 

This theoretical study shows that inclusion of the ex- 
change interaction in the Hartree-Fock approximation to 
the self-consistent solution of the Shrodinger equation 
does allow one to capture the Coulomb blockade effect 
of electron transport. There is no need in the correlation 
effects, at least, for qualitative agreement. The model is 
implemented within the Green's function formalism that 
is well suited to the description of structures of arbitrary 



geometries 
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It was previously tested (without 



the Fock term) on different low-dimensional structures 
and the results obtained were in good agreement with ex- 
perimental data and other theoretical models i^ i^^'^^ An 
attempt to account for exchange interaction was previ- 
ously done using DFT in the local density approximation. 
The conductance however was found to be untrustwor- 
thy: It did not reproduce the 0.7 anomaly measured in 
almost every experiment on quantum point contacts.^— 
This is because the local density approximation of DFT 
lacks the derivative discontinuity and cannot fundamen- 
tally address the transport regime when quantization of 
electron charge occurs-iSiiSrJ^ Another approach that ac- 
counts for the exchange interaction exactly is the Hartree- 
Fock one. It relies essentially on the density matrix that 
takes into account the phase relations between different 
electron states. Accounting for the exchange interaction 
by the density matrix brings qualitatively new physics in 
comparison to the Hartree approach. In the regime of 
a single channel open for transmission, the resonant lev- 
els in few-electron quantum dots become ordered at the 
Fermi energy and conductance oscillations narrow. Both 
of these phenomena are indicative of suppressed screen- 
ing. As potential barriers are imposed between the dot 
and leads, the number of electrons in the dot becomes 
quantized and the conductance shows well defined peaks. 
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FIG. 1: External confinement potential Vext of a quantum 
dot. Only potential below the Fermi energy Ef = is shown. 
The parabolic circular potential confines electrons in the cav- 
ity attached to two semi-infinite leads of width w = 80 nm 
and depth Visad ~ —7 meV. The dot diameter is kept con- 
stant, 320 nm, while the number of electron inside the dot 
is varied. This is accomplished by changing the saddle point 
potential V3. 



To identify the transport regime a different model based 
on the orthodox theory of Coulomb blockade^ was im- 
plemented and calculations for similar but isolated dot 
were performed. The occupancies and conductances for 
both models showed good qualitative agreement, which 
points to the Hartree-Fock approach being adequate for 
describing the Coulomb blockade physics. 

The paper is organized as follows. In Sec. II, the 
Hartree-Fock approach combined with the Green's func- 
tion formalism is formulated. Sec. Ill presents the 
master equation approach to electron tunneling in CB 
regime. The results of numerical calculations and com- 
parison between two models are given in Sec. IV. Dis- 
cussion of obtained results is given in Sec. V followed by 
conclusion. 



II. MODEL 

Consider a two-dimensional open quantum dot at- 
tached to semi-infinite leads. Such a system is routinely 
fabricated in the two-dimensional electron gas by the 
split-gate technique. The gates induce an external elec- 
trostatic confinement Vext that is well described by a har- 
monic potential, Fig. [TJ 

The Green's function of the system in equilibrium is 
defined aai 

(1) 

where m* is the electron effective mass. The second term 
is the electron kinetic energy, while last three terms de- 
scribe effects of the leads^ and electron-electron inter- 
action in the Hartree-Fock approximation. The Hartree 



term accounts for classical electron repulsion 

^Hir) - / dr'n(r') ( - ^ = | , 

' AireoerJ ^ \\r-r'\ ^|r-r'|2+4&2 J 

where n(r) is the electron density screened by the mirror 
charges at the distance of 6 = 30 nm from the surface, 
Sr is the dielectric constant, and the integration is per- 
formed over the whole device area including the semi- 
infinite leads. The Fock self-energy can be written in a 
similar fashion 

SF(r,r') = —n(r,r') I , ^ , ^ - | , 

^' ' iTTEnSr ^' ^\^|r-r'| ^|r-r'|2+4&2 J ' 

(3) 

though its meaning is different: it is the exchange poten- 
tial that respects quantum mechanical principle of an- 
tisymmetry of the wave function. One can show that 
the Fock term completely compensate the Hartree one 
for a single electron thus eliminating self-repulsion, see 
Appendix. In order to remove infinite dimension Sj? is 
adiabatically damped to zero in the leads. Thus, the 
dot is treated in the Hartree-Fock approximation but the 
leads are within the Hartree one. Note that spin degree 
of freedom is out of scope of present study. The phase 
relationship among the different states is established by 
the density matrix 

n(r,r') = -^3 J dE {r,r' , E) fpoiE - Ep) (4) 

whose diagonal elements are simply the electron density 
n(r), fpoiE — Ep) is the Fermi-Dirac distribution func- 
tion and Ep is the Fermi energy. Eqs. ©-(HI) are solved 
self-consistently, see Ref. [13 for details. 

The linear conductance through the quantum dot is 
given by the Landauer formula-"^ 

r [ ,^^,^. dfpD (E-Ep) 

T{E) = Tr [rL{E)g''{E)rn{E)g^{E)] (6) 

where T{E) is the total transmission coefficient, Q°'{E) — 
[Q^'iE)]^ and Tl/r{E) = 23 [Y.L/p{E)\ describes cou- 
pling the scattering region with the leads. 

III. MASTER EQUATION APPROACH 

When the coupling between the dot and leads is weak 
electron transport is governed by Coulomb blockade (CB) 
that in turn is well described by the orthodox theory.^'— 
It assumes incoherent sequential tunneling of electrons 
between the dot and electrodes. The approach is gener- 
ally valid if G ^ ^ . Having small coupling between the 
dot and leads allows one to consider the dot and leads 
separately and treat the coupling as a perturbation. The 
Hamiltonian of the system reads 

H = Hdot + H^lJr + Hp/R^ C^) 
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where H, 



dot 



V,, 



Vh is the Hamiltonian 



2m' ^ 

of the isolated dot and H^^J^ is the Hamiltonian of the 
isolated left and right leads. The latter is equivalent to 
lead self-energies in ([T|). The last term in Eq. ([7]) is the 
standard tunneling Hamiltoniar^ (L/R superscripts for 
left/right lead omitted here and below otherwise speci- 
fied) 



(8) 



with the tunneling matrix elements Tik providing the cor- 
responding tunneling rates 
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(9) 



The total rates of electron tunneling to/from the dot 
are the sum of all partial tunneling rates to/from specific 
energy levels Ek 

T+iN) = Y,^kiE)fFDiE-EF)[l~gN{Ek)]im) 

k 

= Y.^k{E)[l- fFD{E - EF)]gN{Ek)ill) 

k 

where gniEk) is the single particle distribution function 
associated with the Gibbs distribution F{Eki, ■ ■ ■ ,Ekj^) 
of N electrons in the dot, 



gwiEk 



F{Ek„ 



ki. 



Z = 



E 



exp 



k 



^^NN 13) 



Here Z is a normalization constant, and fiN denotes the 
chemical potential in the dot occupied by N electrons. 
To calculate tunneling rates F^, the self-consistent solu- 
tion for Hdot^ = E'^ should be found first. Because we 
are not interested in CB details but rather need to test 
predictions of the Hartree-Fock model established in pre- 
vious section the electron-electron interaction in the dot 
is restricted to the Hartree approximation; Vh is given 
by Eq. ^ with integral restricted by the dot area. 

The electron transport in the Coulomb blockade 
regime is described by a kinetic "master" equation for 
p{N), the probability that there are N electrons in the 
doti^ 

F+(iV - l)p{N - 1) + F"(iV + l)p{N + 1) = (14) 
[T-{N)+T+{N)]piN). 

The average current for any of the tunneling junctions is 
calculated as (with L/R superscripts explicitly written) 

Il/r = eJ2p{N) [F+/^(iV) - F-/^(7V)] . (15) 

N 



This readily allows one to obtain the conductance G 
IlV = -IrV. 



IV. RESULTS 

In order to test predictions of the Hartree-Fock ap- 
proach we compare them with the results obtained us- 
ing the Hartree and master equation approaches. The 
Hartree approximation is known to be correct for open 
systems where several propagating states enter and leave 
the scattering region, i.e. G > The master equa- 

tion is, in opposite, well suitable for description of closed 
systems where CB physics dominates, i.e. G ^ ^-^^ 
As a testing platform we consider an open quantum dot 
with one propagating mode available for electrons to en- 
ter and leave the dot. Fig. [U and then will impose the 
tunneling barriers at the dot openings. The barriers are 
accomplished by a potential raise and equivalent to the 
quantum point contact constrictions formed in experi- 
mental setups. 

As a model system we consider an open quantum dot 
defined by split-gates at GaAs-AlGaAs heterointerface. 
Electron band sturcture in GaAs is well described by 
parabolic dispersion with m* — 0.067m; Sr = 12.9. Note 
however that the model is not material specific and might 
also be applied to other two-dimensional structures, e.g. 
Si inversion layer structures. The parameters of the con- 
finement potential are indicated in the Figure [1] 

The Figs. [2Ja),(e),(i) show the number of electrons N, 
conductance and energy structure as a function of the 
saddle point potential Vg calculated in the Hartree ap- 
proximation. Its remarkable features are linear change 
of N and energy level pinning effect j^iSSiSi The latter is 
pronounced as the resonant levels track the Fermi energy 
within ±2TTkBT interval. Fig. [2{i). Each level represents 
the enhanced density of states (DOS) broadened by cou- 
pling to the extended states in the leads. By minimizing 
the total energy those peaks shift causing metallic-like 
behavior of electrons in the dot. The screening ability 
of the dot is greatly enhanced allowing electrons to eas- 
ily change their location in respect to external electric 
field. One may loosely estimate the energy level sepa- 
ration as A = 0.5 meV, Fig. [D^i). The conductance 
through the dot depends on coupling of particular reso- 
nant states with the extended states in the leadsi^i^l 
Because other states in the dot couple less effectively 
they sweep through the Fermi energy causing no effect 
of the conductance but contributing to the occupancy. 
This explains the period broadening of the conductance 
oscillations in the Hartree approximation^^ 

If the Fock term is taken into account the conduc- 
tance and energy structure change dramatically. Figs. 
mb),(f),(j). The energy structure becomes more compli- 
cated with levels near Ef ordered and less pinned. The 
level separation A about doubles in comparison to the 
Hartree approach. The conductance oscillations become 
narrower. Each conductance peak corresponds to the res- 
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FIG. 2: Occupancy TV (a)-(d), conductance G (e)-(h) and energy structure (i)-(k) of the dot calculated in the Hartree, Hartree- 
Fock, and master equation approaches. The Hartree and Hartree-Fock results in panels (a), (e), (i) and (c),(g),(k), respectively, 
corresponds to potential barriers of transmittance Ti, ~ 0.2 inserted at entrance and exit of the dot. Insert in (k) shows the 
DOS as a function of energy that identify resonant level position. Arrows in (a) and (c) mark Vs potential for charge density 
plots in Fig. [S] Temperature 0.5 K. 




FIG. 3: Representative charge densities in the dot calculated 
in the Hartree and Hartree-Fock approaches for similar dot's 
occupancy, see arrows in Fig. 2(a),(c). 



onant level crossing the Fermi energy. The distance be- 
tween successive peaks equals roughly to the potential Vs 
needed to change dot occupancy by one electron. How- 
ever, N still follows linear dependence on Vs- It becomes 
step-like if the potential barriers separating the dot from 
the leads are imposed. The transmission of each bar- 
rier equals to Tf, ~ 0.2 in Vs range of Figs. [2jc),(g),(k). 
The conductance for a weakly coupled dot generally de- 
creases and well separated peaks become clearly visible. 
Inspection of the electron density reveals strong Freidel- 
like oscillations caused by the exchange interaction in the 
Hartre-Fock approach, see Fig. [31 



The Figs. [2lld),(g) show the results of the master equa- 
tion approach for the dot of the same geometry as studied 
above. The dot however is isolated from the leads and 
treated separately. The agreement between master equa- 
tion and Hartree-Fock approximation for weakly coupled 
dot is remarkable: the electron number is quantized and 
conductance shows periodic peaks. Each peak occurs 
when occupancy changes by one, i.e. in between of quan- 
tized N plateau. All these features are manifestation of 
the CB physics. 

V. DISCUSSION 

Incorporation of the electron-electron interaction in 
the Hartree approximation proved to be sufficient if cou- 
pling between scattering region and leads is strong , ^i^'^i^^ 
G > ^ . It was shown that Hartree and DFT approaches 
provide qualitatively similar description of the electron 
transport in open quantum dots.^" Using the fact that 
DFT predicts spin polarization in a quantum wire simi- 
larly to the Hartree-Fock approximation ™ we may con- 
clude that all these approaches are qualitatively equiva- 
lent for description of electron transport in the strong 
coupling regime G ^ However, in the opposite 

regime of weak coupling, the Hartree approximation is 
not accurate. That regime is governed by Coulomb block- 
ade physics when the electron density is quantized and 
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therefore the Hartree-Fock approach should be rather 
employed. Incorporating the density matrix into the the- 
ory is essential because no features of charge quantization 
and CB are found otherwise. Even though the exchange 
interactions might be treated within DFT in the local 
density approximation or Slater approximation^^ they 
fail to address charge quantization and CB physics. 

The important feature of the present model is absence 
of any adjustable parameters. In this respect, it might 
be considered as a first-principle approach. The mod- 
els existing in literature use adjustable charging parame- 
ters and/or disregard degrees of freedom of the scattering 
region^"— e.g. the dot is treated as a zero-dimensional 
regioni^ Such simplifications make it difficult to relate 
the employed parameter set with physical processes they 
represent and, as a result, the reliability of the results 
obtained is questionable. Note that some studies con- 
sider the device geometry explicitly but still rely on some 
assumption heavily affecting the obtained results. For 
example, slight modification of the Landauer formula al- 
lowed authors of Ref. to address CB even though am- 
plitude of successive conductance peaks fell off exponen- 
tially, which disagree with experimental fundings. 

While the Hartree-Fock model appears to be powerful 
for the description of electron transport it is still a mean- 
field approximation to an interacting many-electron sys- 
tem. The model misses correlation effects that stem from 
the fact that true ground state wave function is not a 
single Slater determinant Electron correlation interac- 
tions might play an important role as it is shown for the 
Kondo effect and others. ^'^■^^ It worth noting also that 
the Hartree-Fock model is known to underestimate elec- 
tron screeningj2S A possible improvement might be ac- 
complished by treating electrons participating in trans- 
port by some exact method while using the mean-field 
approach for others. 

Numerical calculations in the Hartree-Fock approach 
demand much larger computer resources in comparison to 
the Hartree approximation. The matrix dimensions scale 
as (n X raf' and xm for these approaches, respectively 
(n and m are the numbers of sites in transverse and longi- 
tudinal directions, respectively). In practice, the trans- 
verse space coordinate is transformed into momentum 
space where several basis functions are usedi^ Conver- 
gence of the self-consistent solution in the Hartree-Fock 
approach performs generally worse, especially when lead- 

9 2 

to-dot coupling decreases substantially and G ^ A 
possible improvement might be done as suggested in Ref. 
|23| . where the Fermi-Dirac distribution function in Eq. 
is replaced by the Gibbs distribution. 

Finally, it worth noting that the similar results were 
obtained for the open quantum dots of different geome- 
tries. The numerical calculations were also performed for 
the quantum point contacts, where conductance plateaus 
were found to be somewhat larger if the Hartree-Fock 
approach employed. This signals reduced energy level 
pinning effect in accord with results presented above. 



VI. CONCLUSION 

The Hartree-Fock approximation combined with the 
Green's function formalism represents an unified and 
powerful approach to electron transport that implicitly 
takes the interplay between the statistical and quantum- 
mechanical properties of the confined geometries into ac- 
count. This allows one to address the Coulomb blockade 
physics that dominates electron transport in quantum 
dots weakly coupled to the leads. This is confirmed by 
comparison with the master equation approach for equiv- 
alent closed dots: the number of electrons inside the dot 
is quantized and the conductance shows a peak each time 
the electron number changes by one (disregarding spin). 
The conductance peaks are caused by electron resonant 
transmission when corresponding energy levels cross the 
Fermi energy. 

The present study is limited to the case of spinless elec- 
trons. Work is in progress to include the effect of the spin 
in order to revisit magnetic bound state formation in the 
quantum point contacts. It is also interesting to examine 
multiple periodicity in the Aharonov-Bohm interferome- 
ters as well as the effect of mesoscopic CB observed in 
the recent experiment^ on the open quantum dot. 
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Appendix A: Appendix: Self-interaction correction 
in the Hartree-Fock approach 

In order to understand why the Hartree-Fock approach 
eliminates self-interaction error, let us consider a system 
of N electrons each described by the Hamiltonian hi. The 
ground state energy reads 

N N N 

E:^^{i \h,\ i)+^Yl ^ K^^-j'T') > (Al) 

where J represents Coulomb interaction between elec- 
trons in state i and j 

and K represents the effect of antisymmetrization, i.e. 
exchange, 

J J |ri -1-21 

(A3) 

Since Jia;ia' ~ Kic-ia-' clcctron in state i does not self 
interact. This holds true regardless the form of the wave 
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function. Note that in DFT, the Coulomb term J is the proximate. So cancelation of self-interaction in Coulomb 
same as in the Hartree-Fock but the exchange term is ap- and exchange terms is incomplete. 
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